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Abstract 

We study minimizers of the Lawrence-Doniach energy, which describes equilibrium 
states of superconductors with layered structure, assuming Floquet-periodic boundary 
conditions. Specifically, we consider the effect of a constant magnetic field applied 
obliquely to the superconducting planes in the limit as both the layer spacing s — > 
and the Ginzburg-Landau parameter re = e _1 — > oo, under the hypotheses that s = e a 
with < a < 1. By deriving sharp matching upper and lower bounds on the energy of 
minimizers, we determine the lower critical field and the orientation of the flux lattice, 
to leading order in the parameter e. To leading order, the induced field is characterized 
by a convex minimization problem in M 3 . We observe a "flux lock-in transition", 
in which flux lines are pinned to the horizontal direction for applied fields of small 
inclination, and which is not present in minimizers of the anisotropic Ginzburg-Landau 
model. The energy profile we obtain suggests the presence of "staircase vortices" , which 
have been described qualitatively in the physics literature. 
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1 Introduction 

One characteristic of the high-temperature superconductors is their strong anisotropy, which 
distinguishes them from conventional superconducting materials. Indeed, the cuprate super- 
conductors (such as BSCCO, Ba2Sr 2 Ca n Cu n+ i02n+6+a:, or TBCCO, Tl 2 Ba 2 Ca2Cu30io) are 
crystalline materials with a distinct layered structure. Conduction of superconducting elec- 
trons is favored in cuprate planes, which are separated by insulating layers in the crystal 
structure. Physicists recognize that the most appropriate description of such layered super- 
conducting structures is given by the the Lawrence-Doniach model |LaDoj . which was intro- 
duced in 1971 (well before the advent of the high-Tc materials.) The Lawrence-Doniach en- 
ergy treats the superconducting solid as an array of two-dimensional superconducting planes 
described by a Ginzburg-Landau energy, with discrete coupling of adjacent layers in the 
orthogonal direction. In the limit as the spacing between adjacent layers tends to zero, the 
Lawrence-Doniach energy converges to the anisotropic three-dimensional Ginzburg-Landau 
energy |ChDGtfBaKj . This approximation is adequate for certain high-Tc materials, but not 



for the more anisotropic varieties [Ty] , where it is believed that the layered structure plays a 
strong role. 

Many important questions in superconductivity concern the response of a superconduct- 
ing sample to an applied magnetic field. In particular, superconductors expell an applied 
magnetic field of low strength (the Meissner effect,) and there is a critical value of the applied 
field's magnitude Hqi at which flux lines penetrate the material via vortices. In analogy with 
our previous paper (XBS2J on the three-dimensional anisotropic Ginzburg-Landau model, 
we pose the following question: for a given orientation of the applied magnetic field (with 
respect to the superconducting planes,) what is the value of Hci, and in which direction do 
vortex lines and lines of magnetic flux lie? As in [ABS2] , we choose a triply-periodic setting 
in IR 3 . This represents a superconductor occupying all of space, and periodic vortex lattices 
(as are normally observed in the interior of actual superconducting samples.) 

We now define the analytical setting for our problem. The superconducting planes V n , 
jiGZ, are assumed to be horizontal, with uniform spacing s > 0, 

V n = R 2 x {z = z n := ns}, n G Z. 

On each plane is defined a complex-valued order parameter u n : V n — > C. As in the 
Ginzburg-Landau model, \u n \ — 1 represents the purely superconducting state, and u n = 
the normal state. The superconducting currents interact with a magnetic field, described by 
a vector potential A : R 3 ->■ R 3 via h = V x A. 

We assume periodicity with respect to a fundamental domain fl C M. 3 defined by a fixed 
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basis {v[,V2,V3 = (0,0, L)} with 1 The horizontal period is arbitrary, while the 

vertical period is dictated by the discrete geometry of the system. The geometry of the period 
cell dictates that the distance s between adjacent planes be compatible with the period, so 
we assume 



L 

N 



for N eN. 



We also define the restriction of each plane to Q, 

p n = v n nn = {t x v[ + t 2 v 2 + (o, o, z n ) : t u t 2 g [o, i]}, 



n 



1 N. 



As our analysis will show, in our asymptotic regime the natural period of energy minimizers 
will be much smaller, and the choice of fixed period cell Q will not enter into the geometry 
of the minimizers at all. 

To write the LD energy we introduce some convenient notation. We denote by V = 
(d x , d y ), and A' = (A x (x, y, z), A y (x, y, z)), A' n = (A x (x, y, z n ), A y (x, y, z n )). We write (u n , A) 
as a shorthand for ({u n } n& z, A). Then, the energy in the period domain Q may be written 
as: 

N(s) 



E 



n=l J Pn - 

N 



\(V-zA' n )u n \ 2 + — 2 (K| 2 -l 



dx dy 



n=l 



2A 2 s 2 



iff" A z (x,y,z) dz 



dx dy 



+ 



VxA-L 



dx dy dz. 



Here, e > represents the reciprocal of the Ginzburg-Landau parameter; we will assume 
e < 1, which is typical for type-II superconducting materials. The constant A > represents 
the Josephson penetration depth, and will assumed to be fixed in this paper. (Results on 
the large A limit of the LD minimizers may be found in |ABeB 1 , ABeB2j.) The applied field 
h ex is a given constant vector, which will depend on s and e. 
Next we must define a space of functions for C^ s . 

Definition 1.1. We say {u n ,A) e U if u n G F/ oc (M 2 ;C) for all n e Z, A G iJ/ oc (IR 3 ; R 3 ), 



and there exist functions Uj G H t 



loc\ 



), j = 1, 2, 3 so that: 



Unix + vj) = u n (x)e lLU ^' Zn) , 



holds for all x = (x, y, z) 



i = 1,2, neZ; 
u n+N {x) = u n (x)e^' z "\ neZ; 
A(x + vj) = A(x) + Vujj(x), j = 1,2,3 

= (x',z) G IR 3 . 



(2) 
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That is, the configuration (u n , A) e H is f2-periodic up to gauge transformation. In 
particular, the gauge- invariant quantities such as the energy density and the currents, density, 
and induced magnetic field, 

./ / n T r i—i . ./ i . / \ T f i f Zn A z (x,y,z) dz ~) 

j n {x,y) = Im{u n (V - iA n )u n \, j z {x, y) = -Im Ui n u n - X e Jz ™-i ' > 

p n (x,y) = \u n \, h = V x A, 

are all f2-periodic. 

An important consequence of periodic boundary conditions is the exact quantization of 



the magnetic flux through period planes: see Proposition 2J2 below. Indeed, our analysis of 
minimizers for the Lawrence-Doniach model is dependent on a sharp lower bound derived 
using the flux quantization and a slicing method. 

The high-Tc materials being strongly type-II superconductors, it is natural to consider 
the London limit e — > 0. In addition, the distance between layers being very small, we will 
also consider s — > 0. In addition, we will assume throughout the paper that 

s = e° (3) 

for some a G (0, 1). This choice is interesting for two reasons. First, e s is needed for 
the discrete nature of the Lawrence-Doniach model to be retained in the limit. Indeed, if 
s = 0(e) our result predicts the same leading order behavior for minimizers of Lawrence- 
Doniach as for the anisotropic Ginzburg-Landau model, obtained in |ABS2j . Secondly, our 
results in |ABS3j show that for applied fields imposed orthogonally to the superconducting 
planes, the lower critical field (the smallest field strength at which flux lines penetrate the 
sample) is of order | In e | , whereas for applied fields parallel to the layers, the lower critical 
field is of order |lns|. With the choice (|3]), both critical fields are of the same order of 
magnitude, and so we expect both components of an oblique applied field to contribute to 
the energy at the same scale. Thus, we assume the external applied fields are of the form 

h £ ex = H ex | lne| = a~ l H ex \ lns|, (4) 

with H ex a fixed, constant vector (or 2-form) in IR 3 . In fact, it will be convenient to think of 
vectors H e IR 3 as 2-forms, with the association 

(Hi, H 2l H 3 ) H l( ix 2 dx 3 + E 2 dx z dx x + H z dx l dx 2 . 



Employing a Riemannian metric g = diag (1, 1, A 2 ) to represent the anisotropy in the direc- 
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tions in the Lawrence-Doniach model, we define norms on the 2-form if, 



\\m 9 -i = sup p^ff =\Jm+m+x- 2 m- 

x,Ym 3 l A Is- 1 !-* 
Our main result is the following: 

Theorem 1.2. Assume Q, and ([3]). Lei (u ni£ ,v4 e ) £ H be minimizers of C^ s . Then, 

if* in L 2 (n), 

I l ne l 

■urat/i ii* e 1R 3 constant, and 

lim £ ;^"'"^ e) = J ((1 - «)|#* • e 3 \ + a\\H4g + |/f* - H ex \ 2 ) =: F(H*). (5) 



Moreover, if* minimizes F(H) among vectors 



in 



In Proposition 5.1 , we use Fenchel duality to obtain a more geometrical characterization 
of the minimizing field if*. We show that the optimal if* is in fact the closest point to the 
origin in the closed convex body K + H ex where K is defined by 

fll^-^e3|| fl -i<f, iff/3>^, 
UeK hu + ^Wr-! < * if L/j, < - , 



— a 
2 ' 



That is, the set K has a cylindrical side of height (1 — a), and ellipsoidal caps (see 
figure [TJ) The presence of the cylindrical sides of K provides for an interesting phenomenon 
in the Lawrence-Doniach model which is not observed for minimizers of the anisotropic 
Ginzburg-Landau energy. Indeed, let us consider applied fields H ex with a fixed angle 9 
with respect to the horizontal superconducting planes, 

H ex = (H' ex ,H eXt3 ) = \H ex \ cos 9 el + \H ex \ sin # el, 

allowing the magnitude \H ex \ to vary. 

For sufficiently small 9, | tan#| < A (— -J to be precise, the point if of K + H ex nearest 
to the origin can occur in three different locations, each with its own interpretation. The 
situation is illustrated in figure [2] below. For \H' ex \ = \H ex cos 9\ < if = lies in the set 
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Figure 1: The convex set K. 



K + H ex itself (left part of figure |2j) This represents the Meissner effect, and the magnetic 
field is expelled from the sample (no vortices.) 

^ or 2A|cos6>| < l-^easl < 2A|7in0| ' ^ ne neares t point in K + H ex to the origin lies along the 
vertical edge of the set, and so H is horizontal (center image of figure [2j) Thus, although 
the applied field is not horizontal, flux lines are trapped by the superconducting planes and 
the vertical component is suppressed. This phenomenon has been observed in the layered 
high-Tc superconductors [SKMMWJ, and is known as a "flux lock-in transition." It is an 
example of a physical phenomenon which is not present in the anisotropic Ginzburg-Landau 
model, and used as a justification for the Lawrence-Doniach model in the physics literature. 
We note that recent work by Bauman & Xie |BaX] confirms the flux lock-in phenomenon in 
a different asymptotic regime, as A — > oo with s, e fixed (as in [ABcB2j), in nearly-parallel 
applied fields. 

For \H ex \ > 2 a] line I ' ^ ne neares t point lies on the ellipsoidal lower part of the set K + H ex , 
and so H gradually lifts away from the horizontal direction. (See the right image of figure [2]) 
We note that for larger angles, | tan#| > A (^jr), ^ ne closest point of K + H ex will always 
lie either inside the set or on one of the ellipsoidal caps, and so there is no flux pinning 
phenomenon observed in that case. 

From the form of the set K we may thus obtain a formula for the lower critical field 
H c i(9), the largest value of \H ex \ for which H — is the induced field of minimizers: 

f^rS>, iftan#<A(^), 

jj , q\ I 2Acosy' — V a /' 

n c l\P) = < (l- a ) sind+y/a 2 sin 2 6>-(l-2a)A 2 cos 2 B . f „ . f\- a \ 

2(A 2 cos 2 e+sin 2 0) ' iitantf>A(^— J, 



for 9 G [0, f 
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Figure 2: Determining the induced field for external fields H ex applied at a small fixed angle 
9, with | tan#| < A (^-^J with respect to the horizontal. 

It is natural to ask what shape the vortex lines will have in this regime. In the periodic 
case, and with applied fields on the order | lne|, the result we obtain captures the average 
field over a dense vortex lattice, and so a detailed characterization of the magnetic flux 
lines is not available to us. Moreover, since the order parameters u n are only defined on 
the superconducting planes, the usual description of vortices as zeros of a continuous order 
parameter is not applicable here. Nevertheless, an image of the rough structure of mini- 
mizers may be deduced from our upper and lower bound constructions and the form of the 
leading order term in the energy. First, since s — > 0, to the first approximation the energy 
resembles the anisotropic Ginzburg-Landau energy, except that the relevant length scale is 
s, determined by the coupling term in the energy. Thus, each flux line contributes an energy 
of the order 7r|lns| = a7r|lne| per unit length. The second term in the limiting energy 
F(H) (see ([5]))) arises in this way. However, when a flux line penetrates one of the super- 
conducting planes, it leads to a greater energy cost. Inside each vortex tube of width O(s), 
we observe a two-dimensional Ginzburg-Landau vortex with scale e, yielding an energy of 
7rln(s/e) = 7r(l — a)| lne| per vortex. This is the origin of the first term in F(H). This leads 
us to imagine that the flux lines form a step-like pattern: lying nearly horizontally between 
planes, where the energy per length of each flux line is smaller; then penetrating each plane 
nearly vertically, to produce two-dimensional vortices in each superconducting layer. Thus 
we obtain a first mathematical confirmation of the appearance of "staircase vortices" , which 
were described qualitatively in the physics literature |BuLK[ py] , as illustrated in figure [3] 
below. 

The Lawrence-Doniach model contains several dimensionless parameters, and so there are 
many interesting asymptotic regimes to be considered. Our previous papers |ABSll [ABS3J 
consider applied fields which are parallel to the principal axes, in a simultaneous limit s, e — > 
0. Applied fields parallel to the superconducting planes are particularly interesting, and 
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Figure 3: A staircase vortex: the flux line passes through each plane orthogonally, becoming 
flat in each gap to approximate the limiting direction H*. 

we identify two very different limiting regimes. The first regime, s 2 \H ex \ <C 1, is Ginzburg- 
Landau-like, with a regular array of vortices whose mutual distance 5 ~ \H ex \~ l l 2 s, the ra- 
dius of a vortex core. The second asymptotic regime is s 2 \H ex \ 3> 1, for which the sample ap- 
pears transparent to the penetration of the applied field with no loss of superconductivity, and 
many flux lines pass through each gap. Our hypothesis \H ex \ = 0{\ lne|) = 0( \ lns|) clearly 



belongs to the first case, and the result of Theorem 1.2 is consistent with the Ginzburg- 
Landau-like image of minimizers, although with a more subtle geometry. The leading order 
behavior of minimizers in other asymptotic regimes, both for oblique and parallel applied 
fields, remains open. 

The paper is organized as follows: section 2 contains some preliminary definitions and 
results concerning the variational formulation of the periodic problem, and the use of dif- 
ferential forms. In section 3 we derive an upper bound on the energy of minimizers by 
modifying a basic construction for the anisotropic Ginzburg-Landau energy in [ABS2] . A 
complementary lower bound is obtained in section 4, using flux quantization and a slic- 
ing method. Section 5 is devoted to showing the equivalent geometrical formulation of the 
limiting problem for the normalized induced field, H. 



2 Preliminary results 

As mentioned above, given the anisotropy in the model it will be convenient in most of the 
paper to express various quantities as differential forms rather than vector fields, for ease 
in applying transformations to evaluate integrals. To our vector potential, we associate a 
differential 1-form which we also denote by A, 

A = A(x, y, z) = A x dx + A y dy + A z dz. 
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The induced magnetic field is then a 2-form, h = dA. As above, to express the Lawrence- 
Doniach energy we also define the notation 

d' = d x dx + d y dy, A' n = A' n (x, y) = A x (x, y, z n ) dx + A y (x, y, z n ) dy. 

The energy functional may then be rewritten in an equivalent way as: 



N(s) 



Cl s (u n ,A) = sJ2 I 

71=1 J p n. 

N 



\ \d'u n - iA' n u n \ 2 + -^t(|m„| 2 - l) 2 



n=l 



1 



Pn 



2X 2 s 2 



i r™ A z (x,y,z)dz 

^n—l 



dx dy 



dx dy 



If. |2 

+ - \dA — h e ex \ dxdydz. 
2 Jn 



We continue to denote by (u n , A) G % (as in Definition 1.1) the class of admissible Floquet- 
periodic configurations for C^ s , expressed with A as a 1-form and the differential dojj replac- 
ing the gradient Vujj (j=l,2,3) in ([2]). 

In order to assert the existence of energy minimizers we must fix a gauge in which the 
Lawrence-Doniach energy is coercive. For a constant 2-form, h = | Ylijk e i?fc ^ dx^ A dx k , 
we select as representative, 

i _ 

e ijk hiXjdx k . (6) 



.4 



2 ^ 



That is, if we think of h = (hi, h>2, /13) as a vector field, then A = | h x x. We then say 
(un,A) G "H* if u n G if/ oc (IR 2 ; C) for all n G Z, A G ^^(R 3 ; A 1 (R 3 )), and there exist a 
constant 2-form ft, G A 2 (IR 3 ) such that: 



A = A + A , d*A = 0, A (x + Vj) = A (x), j = 1, 2, 3, x G M 3 ," 



M n (x + u,-) = u n (x)e 



-iA(x,z„)-Vj 



j = 1, 2, £ G 



nGZ 



(7) 



Mn+JV (f) = Un (x)e- iA{£ ^\ x G M 2 , 7i G Z, 



where A is associated to h as in (|6j). Note that (w„,A) G "H* satisfy (|2| with ^(x) 
-Mx)^). 

Lemma 2.1. (a) For any (u n ,A) G % t/iere exists 7 G ff 2 oc (lR 3 ) so that 

(u n e^-^\ A + d^en*. 



Moreover, the constant form h is the componentwise average value of dA in Q. 
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(b) There exists a constant C = C (Q) such that for any (u n ,A) £ 

\\A\\ HH n) < C \\dA\\ L 2 iu) . (8) 
This has been proven in our previous paper [ABS3J. 

An important consequence of Floquet periodic boundary conditions is the quantization 
of flux: 

Proposition 2.2. Let (u n , A) £ "H and assume that for some n, \u n \ > on dP n . Then 

deg(^;dP n )=±- [ dA. 



Pn 



7 lb I o 

Given the estimate ([8]), the Lawrence-Doniach functional is coercive on "H* and the 
existence of minimizers for C^ s in H*, for fixed values of the parameters s = L/N, A, e, h £ ex , 
follows from standard arguments AyS |D]. Minimizers satisfy a system of Euler-Lagrange 



equations which contain singular terms along the planes P n . The regularity of the order 
parameters u n and partial regularity of the magnetic potential A (in the Coulomb gauge) 
has been studied by Bauman & Ko [BaK]. Our treatment will be largely variational, and 
we will only require the following estimate on minimizers, which is proven via the maximum 
principle (see |BaK] ): 

Proposition 2.3. Let (u n ,A) e be critical points of C^ s . Then \u n \ < 1 for all x £ P n , 
for all n £ N. 

3 Upper bound 

In this section we assume that 

s = e a , hl x = H ex \\ne\ = -H ex \\ns\, (9) 

a 

where H ex is a constant 2-form, and < a < 1. We have 

Proposition 3.1. Assuming for any constant 2-form h £ A 2 (M 3 ) and for any s > 0, 

there exists {u s n1 A s ) £ % such that 

C [^t^ ~ 2 1 [(1 ~ a)lhl + aMa + lk ~ Hex ^ 
Remark 3.2. The case h% := h(e~[,e~2) = — which corresponds to the case of a magnetic 
field parallel to the superconducting planes — was treated in Proposition 4-4 °f JABS ffl . thus 
we assume throughout the proof that 

h 3 = h(e-[,e 2 )^0 (10) 
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Step 1, the reference configuration. We recall the construction from Proposition 3.1 of [ABS2J, 
but with s, H ex /a and h/a playing the respective roles of e, h ex , h. We begin by defining a 
reference configuration {v,B) in IR 3 \ (Z 2 x IR). 
Let K = [— |, and / solve 

-A/ = 27t(5 (0 ,o) - 1), inX, 
9 V / = on <9K, 

We let 

B = 2tt x dy, 
so that dB = 2ndx A dy. Then we define v by 

f = e 1 ^, where d(j) = j + B and j = *df, (11) 

where *dx = dy and = —dx. Given p G Z 2 we have /(x) = log |x — p| + f p (x), where 
/ p is smooth in a neighbourhood of p and A/ p = 2n, hence j = d8 p + *df p , where 6 P is the 
angle in polar coordinates centered at p. It follows that j + B = d6 p + u p , where du p = 
and Up is smooth near p. Thus v = e 1 ^ is well defined and for each pGZ 2 we have 

v = e i{8p+ ^ (12) 

where £ p is smooth near p. 

Following the proof of Proposition 3.1 of [ABS2J, we define (v s ,B s ). Let p s be a smooth 
cutoff first defined in the square K equal to 1 outside B(0, 2s) and equal to in £>(0, s), and 
then extended by periodicity to IR 2 . Then extend the definition of p s and v to IR 3 by making 
them ^-invariant and let 

v s = ( Ps v) o $ s , B s = $* S B. (13) 
Here $ s = | lns| 1//2 $o ; where <3?o is linear, independent of s and defined by 

$o(S) = 4 

with (&j)j=i,2,3 a (^-orthogonal basis of IR 3 such that 63 is in the kernel of h — i.e. such that 
h(bs, c) = for all c — and normalized by the requirement -(61, 62) — 2vr. More precisely 

Pi| tf = |6j|fl Vi,j, g(bi,bj) = i^j, 

h(b 3 ,-) = 0, h(b u b 2 ) = 2ira. (14) 



Note that by (10) the vector 6 3 cannot lie in the horizontal plane, hence 

^(eg) i Span^i,^). (15) 
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We have dB s = <$>*(dB) = Q*(27rdx A dy) = \ lns\h/a, and the anisotropic Ginzburg- 
Landau energy of the configuration in a period cell 



C, = $1 



1 1 

2' 2 



following the computations in equations (14)-(16) of [ABS2J, satisfies 
G s (v s , B s ; C s ) 



ICJ I lnsl 2 



< 





h 


+ 


h 






H 


a 


9 


a 


a 





(16) 



where 



G s (v,B;n) 



1 



1 



\(d-iB)v\i + -^(\v\'-iy + -\dB-h 



dV P , 



Note that we integrate all terms with respect to the euclidean volume element dV e . 



□ 



Step 2, Periodicity. The configuration (v s , B s ) can be perturbed to be ^4-periodic: We define 
the set of ^4-periodic configurations H by saying (v,B) G % if (v,B) G Hl oc (WL 3 ; C) x 
A X (]R 3 )) and there exist G i^ 2 c (IR 3 ), j = 1, 2, 3 for which 



5(0 



1,2,3. 



u(- + Vj)= v(-) exp(iujj), i vj 3i 

Then (v s ,B s ) is made ^.-periodic by using rational approximation of the transformation 
matrices. We know that (v s ,B s ) is periodic w.r.t. B s = (6j/| Ins) 1 / 2 ). Then we let B' s be the 
basis such that 

M BiA =[M BsA ] = [\ln s\*M BA ], 

where brackets denote the integer part of each matrix entry, where M AB is the matrix whose 
columns are the coordinates of the vectors of B in the basis A, and where we have set 
B = (bi)i- We define the configuration (u s ,A s ) by 

u s = (p s v) o $ s , A s = $*B, 

where $ s is the linear map which maps the i-th vector in B' s to | lns| 1//2 e*j. It is a perturbation 
of $ s in that $ s = $ s o \l/ s , where maps the i-th vector in B s to the i-th vector in B' s . Its 
matrix in the basis A is 

P s :=M AB ,M BsA = ^| In s\^M AB j ' [{lns^M^ 

and therefore P s tends to the identity matrix J3 as s — > 0. Following the computation in the 
previous step, we define the fundamental period cell associated to the perturbed basis, 
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and conclude that: 

G s {u s ,A s ;Q) G s {u s ,A s ;C s ) G s {v s , B s ; C s ) |lne| 2 



|0| \c,\ 

We note for future reference that 



\C S 



= I lns| 1/2 $ n o * c 



,-) ( a ll"ll.<v 



\h-H ex \ 2 ). (17) 



where $o is fixed and satisfies (15), and \l/ s — > Id as s — > 0. 



□ 



Step 3, Discretization. Next, we define a configuration (u a>n , A s ) for Lawrence-Doniach by 
discretization. We claim that there exists t e [0, s) such that, letting u Sj „(x, y) = u s (x, y, z n + 
t), r t A s (x, y, z) = A s (x, y,z + t),we have 



(19) 



Note that this is not the desired final estimate, as we can only bound Lawrence-Doniach 
with e = s in this way. The full estimate will follow from adjusting u njS in region \u n , s \ ^ 1- 



To prove ( 19 ) we use the following elementary observation: for any continuous, ^4-periodic 



real-valued function /, there exists t 6 [0, s) such that 



fdx dy dz 



Jo n=1 JPn 



x, y, z n + t) dx dy 



dr = s 



N 
n=l J 



x,y,z n + t)dxdy. (20) 



The first equality comes from Fubini's theorem, and the second from the mean value theorem 
for integrals. We apply the above with 

f = \\{d'-iA' s )u s \ 2 + ^{\u s \ 2 -l) 2 , 



obtaining t e [0, s) for which (20) holds. We then define u StTl (x,y) = u s (x,y,z n + t), and 
conclude that 

N f 1 1 

^2 / 2^ d ' ~ iA 's( x,y,Zn + t )) M ^rJ 2 + j^(\u s>n \ 2 - ifdxdy 

n=l 

= I ' \\{d' -ir t A' s )u s \ 2 + ^(K| 2 -Ifdxdy dz. (21) 
For the coupling term, we first calculate, 



'U's.n ^s,n— 1*2 



ij* n , r t A s , z {x,y,z')dz' 

° z n—l 



1 



Z n +t 



Zn-l+t 
Zn+t 



d z (u s (x,y,zyf> A ^y^ dz ')dz 

{d z u s -iA s>z v^^ +t ^A^W d z. 



Zn-l+t 
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Then, applying Cauchy-Schwartz, summing over n, and using periodicity we find, 

N 



n=l 



ij* n , nA s , z {x,y,z')dz' 

z n — 1 



dxdy< \ \d z ii s - iA SjZ u s \ . (22) 



Finally, since A s is ^4-periodic, 



\d{r t A s 



h e I 2 

""ex I 



IdS, - hlJ 2 . 



Putting this together with (21) and (22) we obtain (19) 



□ 



Step 4, Conclusion. We may now prove Proposition |3.1| For simplicity we denote from now 
on r t A s by A s . We modifiy £t s n in the set where |u S) „| 7^ 1- We have 

{\u s>n \ + 1} = P n n $ s _1 (U pez ^(p, 2s) x R) . 



From ( 15 ) and ( 18 ) , the angle between P n and the recipocal image by $ s of e*3 is bounded away 
from 0, and therefore the intersection E PiSjn of each elliptical cylinder $ s _1 (B(p,2s) x M) 
with P n is an ellipse whose aspect ratio depends only on s and is bounded uniformly with 
respect to s. Still from ( 18 ), the diameter of the ellipse is bounded above by Cs\ In s\~ 1 / 2 . For 



small enough s, we may therefore include each ellipse in the disk D™ in the plane P n which 
has radius s and same center. We let T s n = {-D^.plpez 2 (n = 1 . . . , iV) be the collection of 
disks with the same centers but twice the radius. 



Outside these disks, \u S)n \ = 1, and hence (19), (17), (16) allow to bound the energy in 
this region: 



N 

n=l Jp n\T s , 
N 



^\(d' - iA' s (x,y,z n ))u n \ 2 + ^(\u s ,„\ 2 - l) 2 \ dxdy 



1 



„ : l 2 ^ S2 ^Pn\[T s ,„UT s , n _l] 



if z Zn . A s , z (x,y,z')dz' |2 



dx dy 



\ L ]dAs ~ h 



e 12 
ex I 



— £s,s( U s,n, A s ) 

< G s (u s ,A s ) < i (a\\h\\ g + \h- H ex \ 2 ) | lne| 2 + o(\ lne| 



(23) 



To complete the upper bound we modify the functions u Sj7l inside each disk D 2s to incor- 
porate a two-dimensional vortex with core size e, and we call the result of this modification 
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First, we define u s ^ n in w s := D 2s \ D s to be of modulus 1, and such that it interpolates 
between u S)Tl outside D 2s and x — Y (x — c)/\x — c\ on dD s , where c denotes the center of the 
D s . From (11) and the discussion following it, near c = $~ l (p, z) we have, denoting by n 
the projection on the plane xy and X = tt& s (x), 



Weil, ( x I 



X -p 
\X-p\ 



(24) 



hence the gradient of u s ^ n is bounded by C/s in u s , and so is the gradient of x — > (x—c)/\x—c\. 
Since A s is bounded by C\ hxs\ we find that 



/ \(d - iA' s )u s , n \ 2 < —\u a \ < C. 

Ju>. s 



Second, inside D s , we let 



U a , n (x) = p e (\x - c|) 



X — C 



\x — c 



where p e = on [0, e), where p e (t) = (t — e)/e on [e, 2e] and p e = 1 on [2e, +oo]. We have 

/ l\tf-U'.)u., n \'+±(\u. tn \*-l)' (25) 

J D s 



I — 1 



< (1 + |lns| 
<7rln(J) [1 + 0(1)], 



l^)| 2 +^(Pe-l) 2 



i f] i |bl.s|) / ^|A S | 2 



where we have used \x + y\ 2 < (1 + a)|x| 2 + (1 + a 1 )|y| 2 (with a = | lns| x ) in the second 
line. 



It remains to count the number of disks in T s n . From (24) each disk in T s>n contributes 
a degree 1 to the map u S:U . On the other hand, since the configuration (it S) „, A s ) is Floquet- 
periodic with respect to (vi,V2), the total degree in each P n is calculated in terms of the 
magnetic flux, 



deg(w Sjn , dP n 



1 

2^ 



eL4< 



Pn 



In s\ h 
2n a 



\Pn\ \Pn\ 

In e| — —h (e\, e 2 ) = |lne| — —h^. 



2tt 



2tt 



Summing over the planes we obtain, since s Y2 n \Pn\ — 1^1 an d in view of (25) 



N 

£ 

n=l 



J- s.n CP„ 



\(d' -iA' s (x,y } z n ))u s , n \ 2 + -^(\u s , n \ 2 - l) 2 > dxdy 



1 — a 



< N M |lne| 2 (l + o(l)). (26) 
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It remains to estimate the contribution of the excluded ellipses T s n to the coupling term in 
the energy. Since \u sn \ < 1, and from the above the measure of T a n is of order s 2 | lns| (the 
area of an ellipse times their number) and we conclude that 



N 1 
oT 



n=l 



2As 2 



\Uari U. 



i f Zn A z (x,y,z')dz' | o , , x i 

.„_ie Jz »-i y ' \ 2 dxdy < C\ lne|. 



[7s in UT Sin _i] 



Combining with (23) and (26), we obtain the desired conclusion and Proposition 3.1 
proven. 



is 

□ 



4 Lower bound 



Let (u ni A) denote minimizers of C* s for s = 1/N — > 0. From Proposition 3.1, C* s (u n , A) < 
C|lne| 2 , so we immediately deduce an estimate for the induced field h = h s — dA in 
particular 

\\h\\ L 2 < C|lne|. 
In particular, passing to a subsequence we may conclude 

h H inL 2 (fi;A 2 (M 3 )). 



| In e | 

We prove the following lower bound on minimizers: 

Theorem 4.1. Let (u n , A) be minimizers of C^ s , with external field h e ex satisfying Q. Then, 
there exists a sequence N — > oo, s — L/N — > 0, e = s l / a — > and periodic H G Lf oc such 
that 

" H mL 2 (tt;R 3 ), 



I me l 
and 

liminf ^,,7 ' > - 1 - a) H 3 + a H g + H - H ex 2 } 
where H is the componentwise average of H over Q, with components 



Jn 



Hi := f Hi 
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4.1 The planes 

For ui C P n , define 



E n (u n ,A;uj) 



J n (u n ,A;u) 



l\(d'-iA'Ju n \ 2 + -^(\u n \*-iy, 



if* n A z (x,y,z')dz' , 2 



dx dy, 



M n {A-u) 



2s 



Z n -\ Jul 



\dA — h ex \ 2 dxdydz, 



so that 



jV 



< s («n, A) = sJ2[En{u n , A; P n ) + J n {u n , A; P n ) + M n (A; P n )}. 



n=l 



We distinguish those planes for which the contribution to the energy is reasonably small: 
define 

/ e :={ne[l,iV]: E n (u n , A; P n ) + M n (A; P n ) < | lne| 3 } . (27) 

Since 

C|lne| 2 >sJ2 E n(u n ,A;P n ) + M n (A;P n ) > scard (I c t )\ lne| 3 , 



n0 € 



we conclude that 



card (J e c ) < C/{s\lne\). 
We obtain the following version of Theorem 4.1 of |SS] : 

Theorem 4.2. For all sufficiently small e > and for every k £ I e , there exists a finite 
collection B k = B(e) = {-Bfc,i}i=i,...,m fe of closed (three-dimensional) balls such that 

1. Each E>k,i has center a^.i G Pk H and radius r k> i with £\ r k ^ = s. 

2. LetV k = B k nP k . Then, 

{xeP k : |K| - 1| > e?} C Vfc. 

3. Set dk t i = deg(u k ,dB kti r\ P k ), and D k = Y2i=i m k \dk,i\- Then, there exists a constant 
C such that 



D k < C|lne| 2 , and 



(28) 



B k nP k 



2\2 



\(d' - tA' k )u k \ 2 + — (l-\u k \ 2 ) 



+ ^ / \h\ 2 >7rD k \n - (l-o(l)) 



2 ./» L6J 



(29) 
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Proof. This result follows almost directly from the basic two-dimensional vortex ball con- 
struction in Chapter 4 of |SS] . The difference is that the Lawrence-Doniach energy does not 
exactly reduce to the two-dimensional Ginzburg-Landau functional, because the energy due 
to the magnetic field is integrated over the domain Q c R 3 , and there is no control over its 
trace onto the planes P n . Thus, we must modify the construction of |SSj so as to integrate 
h over three-dimensional balls rather than planar disks. 

The only real change in the construction is in Lemma 4.4 of |SS] . in which a lower bound 
is obtained for the energy of an S^-valued map v. Following the notation of [SSJ, we replace 
the definition of X by 

X = I dA, 

Js+ 

with the upper hemisphere of radius r, centered on P n , with exterior unit normal vector 
v. In the proof of Lemma 4.4, this quantity is integrated over the planar disk A r , and by 
Stokes' Theorem the two are in fact equal. Our aim is an estimate of \dA\ 2 on the surface, and 
here there is a minor difference in formula (4.18). By applying Stokes and Cauchy-Schwartz, 
we obtain an estimate of the same form as (4.16), 

^ JdS+=dA r A JS+ r V-L + A/ 

In Proposition 4.3, this estimate is then integrated in r. We define 

F{*,r) = \ [ \(d' - iA')v\ 2 + ± I \dA\\ 

where we recall that A r is a (planar) disk in P n , and is the upper half ball (three- 
dimensional.) The proof of Proposition 4.3 continues as before; the result is an estimate for 
S^-valued maps v on the family of balls {B(t)} teR+ , which replaces (4.15) of |SSj . 

1 [ \(d' -iA')v\ 2 + T -±^ [ |dA| 2 >7r|degMA r )|lii-, 

2 JA r \ui 2 JB+ r 2 

and a similar change in estimate (4.14). The remainder of the construction continues as 
before, except the magnetic energy (denoted by H(A,fl) in |SS] ) will be integrated over 
solid balls, and the coefficient of this term in the final estimate will be linear (and not 
quadratic) in the ball radius r. □ 

We may now state the contribution to the lower bound coming from the vortex balls in 
the superconducting planes: 
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Proposition 4.3. With notation as in Theorem \4 . 4 



liminf^S#%^>i(l-a) ! H- 
e^o |lne| 2 2 V in ' 



(31) 



Proof. We first note that 



E 



|/f <s 2 / |/f < S 2 |lne| 



and so this term in the lower bound (29) may be neglected. 

Using the Floquet periodicity ([7]), we note first that for every k G I e , we have flux 
quantization in each plane given by 



D k > 



i=l 



deg(w fc ; dP k 



271 'p. 



h 



(32) 



We now claim that 

f h= I h 

J{z=t}nn Jp k 

is constant in t, k = 1, . . . , N. Indeed, by the regularity theory for minimizers of the 
Lawrence-Doniach functional (see |BaKj ). h 3 = h(e 1 ,e2) is continuous in Q, and the 2-form 
h is smooth in the regions fl n := Q fl {z G (z n _i,z n )}. Let t G [z n -i, z n ), and define the 
horizontal slab, St '■— Q H {z G (t,z n )}. We then apply the divergence theorem to h in St- 
By periodicity of h, the contribution from the lateral edges of St cancel, and hence, 



h 



h. 



'{z=t}nn J p n 

Applying the above with t = z n _i, we conclude that the flux through each plane is identical, 



and in fact is constant in t, and the claim is verified. In particular, (|32|) implies that 
2nD k > 







i 




1 




> 


h 




fl h > 




[ h 




Jp k 


" Z 


Jo J{ z =t}nn 


" L 


Jn 
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In consequence, 

N-l 



s V E n (u n , A; B n ) + S - [ \h\ 2 > s V E n (u n , A; B n ) + [ \hf 

>s^7rL> fc ln[-j (l-o(l)) 



J- 







s~ 


/ ft 


In 








.e. 



(1-0(1)) 



- (card I e ) In - (1 - o(l)) 

i7 L 6 - 



in 



(i-o(i)) 



□ 



since card/ e = N{1 — o(l)) and ATs = L. We divide by | lne| 2 , and use s = e a to obtain the 
result. 

4.2 A comparison functional 

Take minimizers (u n , A) of C^ s and introduce the interpolation 



Z Z-n—l / \ . Z n Z . 

u n {x,y) H u n _ 1 {x,y)e 



i/* n A*(a;,2/,z)<fe 



e Jz 



A*(a:,y,« / )dz' 



(33) 

We observe that if Wj, j = 1, 2, 3, are the functions associated to (u n , A) <E H (see ([2])) then 
A satisfy the Floquet conditions 



ty(x + Vj) = *(x)e^, A(x + iTj-) = A(x) + j = 1, 2, 3, 



(34) 



with the same Uj. As in |ABSlj we connect the terms of the Lawrence-Doniach energy to 
the covariant derivatives of For convenience we define: 



t = t(z) = — - for z n ^i < z < z n . 



(35) 



Note that t(z) is s-periodic in z. 

The connection between the terms in C^ s and the Ginzburg-Landau energy density of 
\l/ are given in the following lemma. Since the estimates are done in local coordinates and 
there are no metric considerations, it will be convenient to express the quantities as vectors 
rather than forms. 
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Lemma 4.4. Assume (u n ,A) are minimizers of C 6)S , and \1/ is defined as in (33). Then, for 
any point (x, y, z) e IR 3 and for any 7 > we have: 



\(d z -iA z (x,y,z))V(x,y, z)\' 



i f Zn A z (x,y,z)dz 

z n — l 



\(V - iA^\ 2 < (1 + 7 )|^)|(V' - iA'(x,y, z n ))u n \ 2 

+ (1 - t(z))\(V - iA'(x,y, z n ^))u n ^\ 2 ) + (1 + 



|2 



(1 - \V\ 2 ) 2 < 2 (1 - |n n | 2 ) z + (1 - t(z)) (1 - K^ 2 ) 



2\ 2 



(36) 
(37) 

(38) 



+ 



4 



Z 7i — 1 



where g n (x,y, z) are functions with the property 



N 

Y\ / k| 2 <C S 2 |ln S | 



Proof. Let (x,y, z) G f2. The identity (36) follows from an explicit calculation, which is done 
in Lemma 4.2 of [ABSl] . From the same proof we also recall that 

(|*| 2 - l) 2 = (2t(l -t)Re{ulu n ^e ! ^ A ^ z)dz } + t 2 \u n \ 2 + (1 -tf\u n _ x \ 2 - l)' 



-t(l-t) 

< 2t 2 (l -t) 2 
1 



i ff n A z (x,y,z) dz 

^n — 1 



+ t(K| 2 -l) + (l-t)(K_ 1 | 2 -l) 



17 z n — 1 



+ 2[t(|n n | 2 -l) + (l-t)(|« 



n— 1 



< 



iff" , Aafov.aOtte 

^ z n— 1 



+ 2[t(| Mn | 2 -l) 2 + (l-t)(| M „_ 1 | 2 -l) 2 ] 



by convexity in the last line, and using \u n \ < 1, valid for minimizers. Hence (|38j) holds. 
Next, we calculate 

(d x -iA x (x,y,z))t> = (39) 
t(d x -iA x (x,y,z n ))u n + (1 -t)e tIz «-^ Mx ' y ' z)dz (d x - iA x (x, y, z n _i))u n _i j x 

where = # n> i + # n)2 , 

g n ,2 = -H 1 - t)u n -ie / ftydz 
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We then estimate the error terms, 



N 



n=l 



N 



\9n,i\ 2 dz dx dy < ^ 



n=l 
N 



h„ dz' 



dz dx dy 



n=l " n ™ J z 
N 



hy dz' dz dx dy 



^ 2 E 



hy dz' dx dy 



s 2 \\h y \\ 2 L 2 in) < s 2 |lne| 2 = o(l). 



A similar estimate holds for g n2 - Expanding the identity (39), we conclude that, for any 
(x, y,z) G Q and 7 > 0, 

\(d x -iA x (x,y,z))y\ 2 



< (1 + 7) 



t(d x - iA x (x,y,z n ))u n + (1 - t)e^™- lAz(x,v ' z)dz (d x - iA x (x, y, z„-i))u 



71-1 



+ (i + 7~ 1 )kl 2 



< (1 + 7) 



t\(d x - iA x (x,y, z n ))u n \ 2 + (1 - t)\(d x - iA x (x,y, z n _i))u n _] 



+ (l + 7- 1 )|^,| 2 . 

The analogous estimate for (d y — iA y (x,y, z n ))u n is proved in exactly the same way, with 
— h x playing the role of h y . These two bounds yield (37). 



□ 



We now use the interpolated function ^ to bound from below the energy in Q \ U n B n , 



the region exterior to the balls constructed in Theorem 4.2 Since the value of ty(x,y, z) in 
each gap z n _i < z < z n depends on the values u n (x,y) and M n _ 1 (x,y) in the two adjacent 
planes, excising a disk in the plane P n precludes calculating the value of ^ in a cylindrical 
region extending through both adjacent gaps, z n -\ < z < z n+ \. Thus, in our lower bound we 
must consider the smaller domain 

Q:=n\C, (40) 

where 

C := [J{(x,y, z) : 2„_i < z < z n+1 and (x,y,z n ) G B n , }. 



We recall that the anisotropy is represented by the metric tensor g = diag (1, 1, A 2 ), with 

> 

9 ^>Jr 



\u\ 2 = V. k g^ k (uj,Uk) for the 1-form uj = Ujdx^ . 
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As in [ABSlJ we define (for u C M 3 ) the functional 



M,(V,A,u) = / m a (V,A) 



2/3, 



|2\2 




(41) 



with (3 = min{l,A -2 }. We note that M. s is very similar to the classical Ginzburg-Landau 
energy, but the role of the Ginzburg-Landau parameter is played by 2s _1 . Following the 
same steps as Lemma 4.4 of [ABSl] we conclude that 

Proposition 4.5. Let (u n ,A) be minimizers of C^ s , and \l/ as in (33). Then, 



£* s (u n ,A)>M s (y,A;{l) + l [ \h-h ex \ 2 + o(l), 

1 Jn 



and (with Q as in (40) J 



M s (y,A;n)< s y2 / 

n=l Jp 



Pn\B 

N 



\(d' - iA' n ) u n \ 2 + ^(\u n \ 2 - 1) 



n=l J Pn 



1 

2AV 



1 A x (x,y,z)dz 

^n — l 



dx dy 



dx dy + o(l). 



(42) 



(43) 



Proof. Let (x,y,z) be chosen with z n _x < z < z n . From (36) and (37), with 7 = | lne| 4 , we 
have 

1, ,2 1, ,2 
41 1 '19 — 4! ig 

< -(1 + i){t{z)\{V - iA'(x,y, z n ))u n \ 2 + (1 - t{z))\(V - iA'(z,y, ^_ 1 ))^ ri „ 1 | 2 } 



+ 



4A 2 s 2 



if* n A z (x,y,z)dz 

u z n—l 



2 +\(l+l- 1 )\9n\ 2 



< -I t(z)\(V - iA'(x,y, z n ))u n \ 2 + (1 - t(z))\(V - iA'(x,y, 2: n -i)K-i| 



+ 



4A 2 s 2 

N 



«P , A z (z,j/, 2 )dz 



+ Iff 



with g n such that 2~^ n =i Jq | f7n | 2 — ► 0. Dividing (38) by s 2 A 2 , and combining the resulting 
inequality with the above, we obtain: 

iM< + ^(i-l*l 2 ) 2 

< -lt(z)\{V -iA\x,y,z n ))u n \ 2 + (l-t(z))\(V -iA'(x,y,2; n _i)K_i 



4 



+ 



+ 



t(z) (1 - \u n \ 2 ) 2 + (1 - *(*)) (1 - K_!| 2 ) 



2\2 



s 2 A 2 

1 

2A 2 s 2 



i ff n A z (x,y,z) dz 

z n—l 



\9n{ x )\' 
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We may then conclude the pointwise bound, 

1,_ 1, l2 1 



m s (\l/, A) < max 
1 

< - 



|v A *| s , 



I|V|^|| 2 4 ( 



I* 



2\2 



t(z)\(V - iA'(x, y, z n ))u n \ 2 + (1 - t(z))\(V - iA'(x, y, z„_i)K_i 



+ 



+ 



4e 2 



t(z) (1 



\u, 



f + (i - t (z)) (i - \ Un ^\ 2 y 



2X 2 s 2 



i \ Zn AJx,y,z)dz 



+ 'AgJx) 



for all (x, y, z) with z n _\ < z < z n , as long as e < s A / 



The lower bound (42) then follows by integrating over each gap and summing over n, 
using fact that t(z) is s-periodic to combine adjacent terms in the sum and obtain C^ s (u n , A) 



on the right-hand side. For the bound (43) we must be more careful, since we cannot include 



contributions to the Lawrence-Doniach energy from the vortex balls in the planes. Notice 
however that if we integrate only over (x,y,z) 6 Q = Q \ C, the vortex balls in the adjacent 
planes above and below are indeed excluded, and thus integrating and summing over the 
region fl again yields the desired bound (43). □ 



4.3 Lower bound by slicing 



Proposition 4.6. Let be defined as in (33), and M. s as in (41). Then, 



M.(*,A;a) > a m 
s^o |f2||lne| 2 ~ 2 11 " 9 ' 

where the 2- form H is the component-wise average of the limiting field H. 



We derive the lower bound on Ai s (^f, A;Q) using the slicing method of [ABS2J. First, 
denote by C the periodic extension (by integer multiples of Vi,V2,v~3) of the exceptional set 
C to 1R 3 . Then, fix a basis (W%, W2, W3) of 1R 3 consisting of integer linear combinations of 
the basis (v[, V2, V3) which generates the period domain Q. Define the domain spanned by 
this new basis, 

Q ■= {a x Wx + a 2 W 2 + a 3 W 3 : < a d < 1, j = 1, 2, 3}, 



and the perforated domain, 



Q:=Q\C. 



It is then immediately true that: 
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Lemma 4.7. \1/ is Floquet-periodic with respect to the frame (W%, W 2 , W3); that is, there 
exist e #L( r3 )> j = 1, 2, 3 so that 

+ W 3 ) = *(.)e^«, A(- + W^.) = A(-) + ^(.), J = 1, 2, 3 

/ioWs m Moreover, 

j^j^U*, 4 Q) = T^M S (*, A; Q), and (44) 
T^M S (V, A; Q) = T^M S (V, A; Q). (45) 



Note that (44) follows from the fact that the energy density of Ai s is periodic with respect 



to both domains Q and Q. To conclude (45), we write A4 S (^ , A;fl) = J n m s (^f , A)xq, and 
note that (by definition of C) the density m s {^ ,A)xq is also periodic. 

We now foliate the domain Q by planes Ut on which 013 = t, < t < 1. Denote by da g 
the surface measure on the plane Il f in the anisotropic metric g. The (Euclidean) volume 
measure dV = dxdydz on Q is proportional to da g da 3 , with constant of proportionality 
dV(W u W 2 , W 3 )/da g (W l , W 2 ), that is, 

dV = —%—da g da,. (46) 



\Wi A W 2 



fi 



The next step is to estimate the energy in each plane Ilf using the vortex-ball construction. 
In order to do this, we need verify two conditions: first, the planes Il t are perforated by the 
exceptional sets C, and these holes must be accounted for by vortex balls. To do this, 
the total radius of the excluded holes must be sufficiently small. Secondly, the vortex ball 
construction demands that the (two-dimensional) energy in the plane IT t be not too large. 

Lemma 4.8. There exists a measurable subset T C [0, 1] so that: 

1. for every t G T , we have both 

I J {m 2 \d A *\ 2 g + \{1 - l*D (|V|vf|| 2 + - |*| 2 ) 2 ) } < S -§ (47) 



and the total radius of the exceptional set, 

p(U t nC) < s|lns| 4 . (4f 

2. The measure of the complement, \T C \ < C|lns|~ 2 . 
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Proof. First, note that by (44) and Proposition 4.5, we have A4 S (^,A;Q) < C|lns| 2 , so 



by Fubini's Theorem the measure of the set T\ of t G [0,1] for which (47) fails is of order 
s 1 / 2 ! lns| 2 . To estimate the measure of the set T2 of t G [0, 1] for which p(Tlt H C) > s\ lns| 4 , 
we list the cylinders in the exceptional set C = {Ck}k=i,...,K, and note that there are at most 
K — N\ lne| 2 cylinders (one for each vortex ball), and each has radius at most 2s. Then, 
taking p(C k H II t ) = when C k D H t = 0, we have 

/ p(Il t nC)dt = y] [ p(C k r\Tl t )dt<C(2s) 2 N\\ne\ 2 <C's(\ns) 2 . 

Jo k=l Jc k nu t ^ 

The set r 2 of t for which p(H t r\C) > s\ In s| 4 is thus at most of measure | In s| -2 . We conclude 
that T c = 7~i U r 2 has measure \T C \ = 0(\ lns|~ 2 ). □ 

For t G T we apply the vortex ball construction as in |SS] . We have: 

Proposition 4.9. For every t G T , there exists a finite collection of disjoint closed disks 
T> t = {A*} Jg j t lying in the plane U t such that: 



1. the sum of the radii p{V t ) = J2j<=j t p(^j) — I m s \ 10 

2. {x G Tl t ■ ||^| 2 - l| 2 > s 1 / 8 } C V t ; 

3. C n Il t C V t ; 

4. Let d) = deg(^; 9A*.) and D t = £ je j t \d\\. Then, 

1 



r: 



T> t \C L 



do-g > 7rA|lns|(l - o(l)). 



(49) 



The proof of Proposition 4.9 follows from Theorem 4.1 of [SS]. The primary difference in 



the two proofs that the exceptional set C PI IT must be included in the collection the initial 



balls (B in the notation of [SSJ.) For t G T, the total radius of the balls satisfies (48), and so 
including these sets does not affect the remainder of the proof. The density m s (^/,A) may 
be substituted for the usual Ginzburg-Landau energy density since |^| > 1 — s 1 / 8 outside of 
the initial balls (Bo.) 



We note that by conclusion 1 of the Proposition 4.9 and the periodicity of \^/\ we may 



(if necessary) translate the period domain Q in order that no vortex ball in U t intersects 
the boundary. As a consequence, we may define the degree of \1/ along this boundary, 



S. Alama, L. Bronsard, & E. Sandier 



27 



deg( 1 I r ; dU t ). By Floquet periodicity and Stokes Theorem, the degree quantizes the flux of 
the magnetic field h = dA through the surface IT t : 



/ h = 2tt deg(^; dU t ) < 2itD t . 
Ju t 



Thus, the lower bound (49) implies 

1 

m s {V,A) + 

>v t \c 



m s (ty,A) + -r 2 \dA\\ 



da g > l -yj^ h)\lns\(l-o(l)). 



We now integrate the lower bound (|51|) over t £ T, using the identity (|46|): 

1 



m s ($, A) + -r \dA 



> 



> 



t£T JV t \C 

\Q\ 



dV 

[m s ^,A)+r 2 \dA\ 2 g ] 



\Q\ 




2|WiAW 2 | \J T Jn 

i \Q\ 

2 \Wi A W 2 \ 



W 1 AW 2 
h ) |lns|(l-o(l)) 



-da g dt 



I [ h) |lns|(l-o(l)), 
Jo Jn t J 



since by Lemma 4.8 we have 





< \T C \ 







Now, since we have 
• i 



1 



h= [ I [ h(W 1 ,W 2 )da 1 da 2 dt = — / h(]\\.]\' 2 )dV. 
Jilt Jo Jo Jo \Q\ 



the right-hand side simplifies to 



1 f h{W!,W 2 ) 



dV. 



2 J Q \Wi AW 2 \ g 

Finally, we remark that the term r 2 |cL4| 2 above may be neglected. Indeed, 

/ r 2 \dA\ 2 g < Cr 2 [C^(u n , A) + \Q\\h ex \ 2 ] < o(l). 

J Q 

Therefore, we rewrite the lower bound as: 

M S {V,A;Q)>1- [ ,yV^7 dV | In a| (1 - o(l)). 



2J Q \Wi A W 2 \ g 



(50) 



(51) 
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Writing h as the constant 2-form representing component-wise average of the 2-form h, we 
— > H, the component-wise average of the limiting field. Note that by periodicity, 



have 



|lne| 



the averages are the same for any period domain, and so (by Lemma 4.7) we conclude 



, Ms(*,A;Q) , Ms(*,A;Q) a H(W u W 2 ) a H(W u W 2 ) 

|n||lnse| 2 ™ |Q||lne| 2 ~ 2 \W X A W 2 \ g ~ 2 |WiWW 2 |p ' 

Taking the supremum over all such bases (W 7 !,^,^) which are integer multiples of the 



periods we obtain the desired lower bound, and Proposition 4.6 is complete. 



4.4 Finishing the proof of Theorem 1.2 



Putting together the results above, and using s = e a , we have: 



CUu n ,A;tt) 



liminf 

lit lneP 



> F(H) = - [(1 - a)|/f • e 3 | + a||^|| ff + \H - H ex \ 2 ] . 



Applying the upper bound from Proposition 3.1 with h = H* the minimizer of F(H) over 
ffeR 3 produces the complementary inequality, 

CUu n ,A; fi) 



lim sup 



s— >0 
e->0 



Ifilllnel 2 



< F(^) < F(H). 



In particular, = H = is constant, and the normalized energy converges, proving 

(|5]). This completes the proof of Theorem 1.2 



5 Minimizers of the limit energy 

In this section we use convex duality to derive a more geometrical characterization of the 
limiting normalized induced field, H. In this setting, it will be convenient to represent the 
2-form as a vector field, with the appropriate norms as defined in the Introduction. 



Proposition 5.1. Let 

F(H) 

Then the Fenchel dual is 

F*(H) -- 



1 2 1 — a a 

~\H — H prr \ + \H-t \ H — \\H \\ n 

2 1 1 2 1 1 2 " 119 



l\H\*-hH, 12 



-oo, 



■■) i 11 1 .v ■ if H II, , G l\ . 

if H — H ex £ K, 
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where K is defined by 



l-g 
2 ' 



\U - ^esWg-i < f, ifU 3 > 
UeK ^ \ wu+^e^ < a j/^ 3 <_(!=«), (52) 

where we denote by U' = (Ui, U2, 0) and U 3 = U ■ e 3 . Thus, if H* G IR 3 minimizes F , then it 
also minimizes F* , and F(H*) = —F*(H*). 

Proof. Let U = H — H ex , equivalently, we minimize 

G(U) = ^\U\ 2 + V(U), V(U) = ^-^-\U 3 + H ex , 3 \ + |||C/ + H ex \\ g . 
By convex duality, 

min F(H) = min G(U) = - min -\U\ 2 + V*(-U), 



um 3 um 3 2 

where \&* is the Fenchel dual of 

qf*(U) = sup [U-V-%(V)\ 



- inf [-U -V + V(V)] 

veR 3 



—U ■ H ex — inf 



1 — a , , a ,, 
_ f/ . y+ _|K.e 3 | + -||y|| 9 



If there exists V G IR 3 with -£7 • V + ■ e 3 \ + § \\V\\ g < 0, by homogeneity it follows 

that the infimum is —00, and then ^>*(U) = +00. On the other hand, if U G K, 

K := |[7 G M 3 : -17 • V + ^^1^ • e 3 | + > 0, for every 7 G M 3 ,| (53) 

then the infimum is attained with V = 0, and \?*(Z7) = —U ■ H ex . 

We claim that K defined above coincides with (52). We divide into the three cases as 
above. 

Case U 3 > First assume U G K. Let Vq G IR 3 be the vector with ||Vo|| 9 = 1 which 

attains the supremum, 

1 — a , 1 — a , , 1 — a x 

\\U - -z-e 8 || a -i = sup (U - -— e 3 ) ■ V = (U - -— e 3 ) ■ V . 

1 ni s =i 1 * 

Since U 3 — > 0, clearly the sup is attained with V03 = Vo • e 3 > 0. Thus, 

1 — OL 1 — OL 1 — OL OL 

W ~ — e 3 \\ g -i = (U- -^e 3 ) . Vo = U-V - -^\V 03 \ < |, 
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as desired. 

Conversely, assume \\U — ^pe 3 || s -i < |. Then, for any V G M. 3 , 

ol ,, ,, , 1 — a N 1 — a,, 

-\\V\\ 9 >(U- — e 3 ) • V > 17 • V - —\Vs\, 

that is, U e K. 

Case ?7 3 < — (— ^p): As above, assume U G K, but now let Vq, \\Vo\\ g = 1, attain the 



supremum m 



1 — a . 1 — a N . 1 — a N 

1 17 + —z-eaWg-i = sup (C/ + — — e 3 ) • V = (U + ~^—e 3 ) ■ V . 
1 ||V|| fl =l ^ ^ 



Now we must have V03 = Vq ■ e 3 < 0, and hence 

\\U + -g— e 3 || fl -i = (C/ + -g— e 3 ) • Vo = U ■ V - — ^ — | V 03 1 < ^ 
as desired. The converse is proven exactly as in the previous case. 

Case - (±=2) <U 3 < (i=«): First assume U <E K, and let V = K) = (Hi, V02, 0) with 
||Vo||g = 1 attain the supremum in 

\\U'\\ g -i= sup U' ■ V = U' ■ Vq. 

v'em 2 x{o} 

\\V'\\g=l 

Then, 

WU'L-i = U' ■ V ' = U ■ V - ^\V 03 \ < 



2 



as desired. 

Conversely, if llC/'IL-i < ~, for any V G IR 3 we have 



is — 2 
a 



-\\V\\ g >U'-V = U-V-U 3 V 3 

>U - V — \u 3 \ \v 3 \ 
> u ■ V - ^-^\v 3 \. 



This concludes the proof of the claim (52). 
We conclude that 

{+00, ifU^K, 
and reverting to the original variable H we conclude that 



F*{H) 



\\H\ 2 — \\H ex \ 2 , if H — H ex G K, 
+00, if H - H ex 4 K. 



□ 
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